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Abstract. We study non-equilibrium statistical mechanics of a Gaussian dynamical system 
and compute in closed form the large deviation functionals describing the fluctuations of the 
entropy production observable with respect to the reference state and the non-equilibrium 
steady state. The entropy production observable of this model is an unbounded function on the 
phase space, and its large deviation functionals have a surprisingly rich structure. We explore 
this structure in some detail. 
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1 Introduction 

In this paper, we prove and elaborate the results announced in Section 9 of [JPR]. We consider a dynamical 
system described by a real separable Hilbert space K, and the equation of motion 

^xt=Cxt, xo & 1C, (1) 

at 

where £ is a bounded linear operator on K,. Let £> be a strictly positive bounded symmetric operator on K, 
and (X, (jJd) the Gaussian random held over 1C with zero mean value and covariance D. Eq. (1) induces a 
flow (j)c = {(l>c} °tt ^tttl our starting point is the dynamical system (X, (its detailed construction 

is given in Section 2.1). We compute in closed form and under minimal regularity assumptions the non¬ 
equilibrium characteristics of this model by exploiting its Gaussian nature. In particular, we discuss the 
existence of a non-equilibrium steady state (NESS), compute the steady state entropy production, and study 
the large deviations of the entropy production observable w.r.t. both the reference state ujd and the NESS. 
To emphasize the minimal mathematical structure behind the results, in the main body of the paper we 
have adopted an abstract axiomatic presentation. In Section 3, the results are illustrated on the example 
of the one-dimensional harmonic crystal. Eor additional information and a pedagogical introduction to the 
theory of entropic fluctuations in classical non-equilibrium statistical mechanics, we refer the reader to the 
reviews [RM, JPR]. 

There are very few models for which the large deviation functionals of the entropy production observable 
can be computed in a closed form, and we hope that our results may serve as a guide for future studies. 
In addition, an important characteristic of a Gaussian dynamical system is that its entropy production 
observable is an unbounded function on the phase space. This unboundedness has dramatic effects on the 
form and regularity properties of the large deviation functionals that require modifications of the celebrated 
fluctuation relations [ECM, ES, GCl, GC2]. Although this topic has received a considerable attention 
in the physics literature [BaCo, BGGZ, BJMS, Fa, HRS, Vil, Vi2, ZC], to the best of our knowledge, it 
has not been studied in the mathematically rigorous literature on the subject. Thus, another goal of this 
paper is to initiate a research program dealing with mathematical theory of extended fluctuation relations in 
non-equilibrium statistical mechanics, which emerge when some of the usual regularity assumptions (such 
as compactness of the phase space, boundedness of the entropy production observable, smoothness of the 
time reversal map) are not satisfied. 

The paper is organized as follows. In Section 2.1 we introduce Gaussian dynamical systems. In Section 2.2 
we define the entropy production observable and describe its basic properties. In Section 2.3 we introduce 
the NESS. Our main results are stated in Sections 2.4 and 2.5. The entropy production observable is defined 
as the phase space contraction rate of the reference measure ujd under the flow (pc, and in Section 2.6 we 
examine the effects of a perturbation of the reference measure on the large deviation theory. In Section 3 
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we illustrate our results on two classes of examples, toy models and harmonic chains. The proofs are given 
in Section 4. 

The focus of this paper is the mathematics of the large deviation theory of the entropy production observ¬ 
able. The physical implications of our results will be discussed in the continuation of this paper [JPS]. 
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2 The model and results 

2.1 Gaussian dynamical systems 

In order to setup our notation, we start with some basic facts about classical Gaussian dynamical systems. 
We refer the reader to [CFS] for a more detailed introduction to this subject. 

Let r be a countably infinite set and 


X = {x = (x„)„gr I x„ e K} = 


Eor X G X and J C T, we denote xj = {xi)i^i G K^. Let I = ((n)raGr be a given sequence of strictly 
positive numbers such that ^ (we shall call such a sequence a weight). Then 



is a metric on X and (X, d) is a complete separable metric space. Its Borel cr-algebra F is generated by the 
set of all cylinders 


Ci{B) = {xGX\xiGB}, 


where / C T is finite and C is a Borel set. 

Let V and w be two Borel probability measures on X. We shall write <C a; when v is absolutely continuous 
w.r.t. w. The corresponding Radon-Nikodym derivative is denoted by 



We will also use the notation 




The two measures v and w are called equivalent, denoted ~ w, if they are mutually absolutely continuous, 
j.e., u! V and v ut. We adopt the shorthand ^(/) = fdv. The relative entropy of v w.r.t. w is 
defined as 


Ent(i/|w) 



if ^ <C w, 


otherwise. 


( 2 ) 


Throughout the paper we adopt the convention log x = —oo for a: < 0. 
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We recall that Ent(^|a;) < 0, with equality iff ly = uj. For a G M, the relative Renyi a-entropy of ly w.r.t. oj 
is defined as 

f logo; if^<a;, 

EntQ(z^|a;) = < 


otherwise. 


We denote hy /C C X the real Hilbert space with inner product 

{x,y) = '^Xnyn (3) 

Tier 

(JC = £g(r)), and by {(5„}„gr its standard basis. The matrix elements of a linear operator A on ^^(r) 
w.r.t. this basis are denoted by Anm = (i5„, A5m)- 

Let X/, Xj* C X be real Hilbert spaces with respective inner products 

( 2 ;, y)i = X! ^nXnyn, {x, y)i. = ^ 

nGr nGr 


(Xj* is the dual of X; w.r.t. the duality (3)). Clearly, 

cKcXiCX, 


with continuous and dense inclusions. All the measures on {X,T) we will consider here will be concen¬ 
trated on X/. 

Let D he a bounded, strictly positive operator on 1C. The centered Gaussian measure of covariance D 
on (X, X") is the unique Borel probability measure uj^ specified by its value on cylinders 


ojd{Ci{B)) 


y^det(27rD/) Jb 


3 ^^dj 


where Dj = je/- The measure ujb is also uniquely specified by its characteristic function 


X’^ ^ y x{y) = ( ChjJd{x) = e 

'j 3C 


The bound 


f MfduJD{x)= f ^Z„a;^dw_D(a;) = 

•'^ner ner 


\D\\ 


(4) 


implies that ujd{X \ Xi) = 0, f.e., that lod is concentrated on X/. 

Let T be the real vector space of all trace class operators on JC and ||r|| 1 = tr((T*T)^/^) the trace norm 
on T. The pair (T, || • |ji) is a real Banach space. By the Feldman-Hajek-Shale theorem, two Gaussian 
measures ujdi and on (^i are equivalent iff T = G T. In this case, one has 




(5) 


Ent(a;Dja;£,J = ^ti {DiT{I + DiT) - ilogdet(/ + DiT). 

Note that det (/ + DiT) = det = det{Dl^^> 0. 

Let £ be a bounded linear operator on JC such that £*X^ C X;*. It follows that C has a continuous extension 
to X/ which we also denote by C. For x G X and f G K we set 

{ e*^x if X G X;, 

(6) 

X if X ^ X;. 

The map {t, x) 1 —>■ (jy^n^x) is measurable and (pc = {(pc}*^^ ^ group of automorphisms of the measurable 

space (X, B) describing the time evolution. We shall call pc the dynamics generated by C and (X, Pc,^d) 
a Gaussian dynamical system. Note that for w/j-almost all x G X, Pci^) ~ g*^x for all f G M. 
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2.2 Entropy production observable 

Our starting point is the dynamical system (X, (j), oj), where cj) is the dynamics on X generated by C and uj 
is the centered Gaussian measure with covariance D (from now on, C and D are fixed, and we shall omit 
explicit reference to them). The measure uj is sometimes called the initial or the reference state of the 
system. Observables are measurable functions / : X —> C. They evolve according to 

Mx) = fo4>\x). 

The expectation of an observable / at time f G K is given by 

^t{f) = w(/t) = J Mx)duj(x), 

where w* = w o (j)~* is the centered Gaussian measure on (X, X") with covariance 

Dt = 

Dt is a bounded strictly positive operator on ^^(r) and uJt{Xi) = 1 for all t. By the Feldman-Hajek-Shale 
theorem, the two measures oot and w are equivalent iff Tt := — D~^ G T. We shall assume more: 

(Gl) The map K 9 f i—>■ Tt G T is differentiable at f = 0. 


As will be seen later, this condition implies that the function f i-9 Tt is differentiable for all t. The entropy 
production observable (or phase space contraction rate) for (X, (/>, w) is defined by 


cr{x) =—£^^\^{x) xGiC. 

A simple computation shows that (cf. (37)) 

cr(a;) = {x, ^x) - tr(Z)^), 


(7) 


where 


and the derivative is understood in the sense of T (in particular, G T). Since T is continuously embedded 
in the Banach space of all bounded operators on IC, we have 





Remark. If A is a self-adjoint element of T, then the quadratic form (x, Ax) has a unique extension from K. 
to an element of (X, dw). With a slight abuse of notation, we shall also denote this extension by {x, Ax) 
(see Lemma 4.1 below for a more precise statement). Thus, the entropy production observable (7) is a 
continuous function on JC and an integrable function on X w.r.t. the measure w. 


Proposition 2.1 Suppose that (Gl) holds. Then: 


(1) The function K 9 f i—>■ dt G T^(X, dw) is continuous. 

(2) = /o <X-s ds holds as the Riemann integral of a continuous T^(X, duj)-valued function. It also 
holds for Lo-almost every x G X as the Lebesgue integral of a real-valued function. 

(3) The function M 9 f i-9 g T^(X, dw) is and 


d 

dt 




cr-t- 


(9) 
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(4) = tr((r(Dt — D)) and in particular uj{a) = 0. 

(5) Ent(wt|w) = - J*ujs{<j)ds. 

In specific examples, it may happen that only finitely many matrix elements <;nm are non-zero, and in this 
case the map x i—?> cr(x) is continuous on X. The function a is bounded only in the trivial case a — 0. Note 
that a = 0 iff uJt = uj for all f; this follows, for instance, from the cocycle property (38). 


2.3 Non-equilibrium steady state 

Our next assumptions are: 

(G2) There are some numbers 0 < m < M < oo such that m < Dt < M for all f € M. 
(G3) The following strong limits exist: 

s - lim Dt = D±. 

t—^dzoo 


It is clear that m < D± < M, and CD± + D±C* = 0. In what follows, we set 

(5 ^ 

M — m 

Let uj± be the centered Gaussian measure on (X, D) with covariance D±. 
Proposition 2.2 Suppose that (G1)-(G3) hold. Then: 

(1) For any bounded continuous function / : X —> M, 

lim wtif) =UJ±{f). 


( 10 ) 


(2) cr € L^(X, dw±) and 


u}±{a) = lim wt(cr) = tr(<?(L>± - D)). 

t^±oo 


Note that 


w+((t) = lim - 

t^OO t 


Ws(cr)ds 


lim -Ent(a;t|a;). 

i—>-oo t 


We shall call the NESS and the non-negative number w+(cr) the entropy production of (X, ui). 


2.4 Entropic fluctuations with respect to the reference state 

Time reversal invariance plays an important role in non-equilibrium statistical mechanics, and in particular 
in formulation of the fluctuation relations. Hence, we shall also consider the following hypothesis: 

(G4) There exists a unitary involution "d : K, ^ IC such that 'd{Xi) C X;, dC = —Cd, and 
■dD = Dd. 

This assumption implies that D_t = dDtd for all t G K, and thus D_ = dD^d and w_|_ = a;_ o d. 
Moreover, it follows from Definition (8) that r?? = —This in turn implies that tr(D<j) = 0 and 

cr(x) = uj+{a) =-u}-{a). (11) 

For simplicity of notation and exposition, we shall state and prove our main results under the time reversal 
invariance assumption, which covers the cases of physical interest. With a minor modifications of the 
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statements and the proofs, most of our results hold without this assumption. We leave these generalizations 
to the interested reader. 

The relative Renyi entropy functional, which is defined by 

et(a) = Enta(a;t|a;) = loga;(e“^‘"ti“), (12) 

is a priori finite only for a € [0,1]. To describe its properties, we introduce the sets 

</i = |o;€M|T) ^ + ckTt > 0|, t G. M, 
and denote by C± the open upper/lower half-plane. 

Proposition 2.3 Suppose that (G1)-(G4) hold. Then: 

(1) Jt = {—St, 1 + St) for some St > S and J-t = Jt- 

(2) The function a i— et{a) is finite on the interval Jt and is equal to +oo for a ^ Jt- Moreover, this 
function is convex, extends to an analytic function on the cut plane C+ U C_ U Jt, and satisfies 

ei(0) = et(l) = 0, e;(0)<0, e;(l) > 0. (13) 

In particular, et{a) < Ofor a S [0,1] and et{a) > 0 otherwise. 

(3) The finite time Evans-Searles symmetry et{a) = et{l — a) holds for all t and a. 

We now study the statistical properties of trajectories as f —> -|-c». The intervals Jt do not necessarily form 
a monotone family, and we define the minimal interval 

J = liminf = M Pi Jt. 

t-tao ^ ' < 

T>0 t>T 

Clearly, one has J = {—S, 1 -I- S), where 5 = liminft_>oo St > S. 

Theorem 2.4 Suppose that (G1)-(G4) hold. 

(1) The limit 

e{a) := lim -et{a) (14) 

t—^+oo t 

exists for a € J_. Moreover, the function e{a) is convex on the interval J_ and satisfies the relations 
e(0) = e(l) = 0, e^(0) = — a;+(cr) < 0, e^(l) = u;+(ct) > 0, e(l — a) = e(a). (15) 

(2) The function e{a) extends to an analytic function on the cut plane C+ U C_ U Jj and there is a unique 

signedBorel measure v with support contained in K \ J such that f |r|“^d|z/|(r) < oo and 

e(a) = — [ log fl — —) dt^(r). (16) 

Jr ^ rJ 

(3) The Large Deviation Principle holds in the following form. The function 

I{s) = sup (as — e(—a)) 

-aGJ 

is convex, takes values in [0, oo], vanishes only at s = w_|_ (cr), and satisfies the Evans-Searles symmetry 
relation 

I{—s)=I{s) + s forsGM.. (17) 

Moreover, there is e > 0 such that, for any open set J C (— (cr) — e, a;+ (cr) -f e), we have 

lim - log w (lx G X - [ (Ts{x) ds G = — inf I{s). (18) 

i^oo t t Jq } j 
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(4) The Central Limit Theorem holds. That is, for any Borel set i? C K, w have 

rt 


lim uj 

t—foo 


X G X 




da; 

s/2'Ka' 


where a = e"{!). 

(5) The strong law of large numbers holds. That is, for uj-a.e. x G X, we have 


1 /■* 

lim - / (Ts(a;) ds = a;+(cr). 
t-!-00 t Jq 


(19) 


Remark 1. In general, the two limiting measures uj- and oj+ are distinct. This property is closely related 
to the strict positivity of entropy production. In fact, it follows from the second relation in (11) that if 
oj- = UJ+, then = 0 as well as w_(tT) = 0, while any of these two conditions imply that the 

function e(a) vanishes on [0,1] and, hence, identically in view of analyticity. 

Remark 2. The representation of e{a) as a logarithmic potential of a signed measure is somewhat surpris¬ 
ing, and its mathematical and physical significance remains to be studied in the future. The measure o is 
related to the spectral measure of the operator Q (see the proof of Theorem 2.4 for more details). 


Now let {tn} C K+ be a sequence such that St„ —>■ 6. We define J = {—S, 1 -I- S). Note that, by 
Proposition 2.3 (1), we have <) > (5. In the case when 6 coincides with 6 = limsupj_^go 6t, we write J 
instead of J. 

Theorem 2.5 Suppose that (G1)-(G4) hold and {tn} C K+ is a sequence satisfying the above hypothesis. 
(1) Let Q = d]/^{DZ^ - DZ})d} 1^. Then 

1 1 

( 20 ) 


-= < Q < 

5 “ l + (5 


Furthermore, since the function g{z) = z ^ log(l — z) is analytic in the cut plane C \ [1, oo), the 
operator-valued function 

E{a) = -aD]!'^g{aQ)D]!'^, 


( 21 ) 


is analytic in the cut plane C+ U C_ U J. 
(2) For a G J, the following relation holds: 


e{a) := lim —et„{a) = tT{E{a)<;), 

n—>-oo tn 


and if a GM. is not in the closure of J, then 


1 


lim sup — (a) = cx). 

tn 


( 22 ) 


(23) 


n—foo t'n 


Moreover, the function e{a) is convex on the interval J and satisfies relations (15). 

(3) The Large Deviation Principle holds in the following form. The function 

I{s) = sup (as — e(—a)) 

-ctGJ 


(24) 


is convex, takes values in [0, oo], vanishes only at s = a;+ (cr), and satisfies the Evans-Searles symmetry 
relation (17). Moreover, for any open interval C K, we have 


lim — log OJ 

n—>oo tn 


X G X 


-i 

tn Jo 


(Ts{,x)ds G J \ ) = - in^/(s). 


(25) 
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Remark 1. The functions e{a) constructed in Theorem 2.5 coincide with e{a) on the minimal interval J. 
Moreover, by Part (2) of Theorem 2.5, the functions e constructed for different sequences {t„} must coin¬ 
cide on the common domain of definition. 

Remark 2.lfS = oo, then e(a) = e{a) = 0 for a S K. 

Remark 3. The local Large Deviation Principle described in Part (3) of Theorem 2.4 is an immediate 
consequence of the local G^tner-Ellis theorem (see Appendix A.2 in [JOPP]). The global Large Deviation 
Principle described in Part (3) of Theorem 2.5 cannot be deduced from the Gartner-Ellis theorem. Our 
proof of the LDP exploits heavily the Gaussian structure of the model and is motivated by Exercise 2.3.24 
in [DZ], see also [BLL, BLR, BD] for related results. 


2.5 Entropic fluctuations with respect to the NESS 

We now turn to the statistical properties of the dynamics under the limiting measures uj± . In view of the 
time-reversal invariance (G4), it suffices to study the case of one of these measures, and we shall restrict 
ourselves to uj+. Let us set (cf. Part (2) of Proposition 2.1) 

et+(a) = logw+(e““^“‘i‘") = log w+(e““-^o= loga;+(e““-^o 

where the last relation follows from the invariance of a;_|_ under the flow (/)*. Note that, a priori, et+{a) 
might not be finite for any a ^ 0. 

Theorem 2.6 Suppose that (G1)-(G4) hold. Then: 

(1) For any f € K, the function K 9 a i—>■ et+{a) G (—oo, -foo] is convex. 

(2) The set 

Jt+ = {a G M I - aTt > O} (26) 

is an open interval containing (—6,6), and the function et+{a) is real analytic on and takes 
value -foo on its complement. 

(3) Let be the interior of the set 

T>0 t>T 

Then J'*’ is an open interval containing {—6, 6). Moreover, for a G , the limit 

e-t(«) = lim Tet+{a) (27) 

)-oo t 

exists and defines a real-analytic function on jf'. Finally, if a is not in the closure of , then 

limsup -et+{a) = -foo. (28) 

t—>oo t 


(4) The Large Deviation Principle holds in the following form. The function 

I^{s)= sup (as —e+(—a)) 

-aG 

is convex, takes values in [0, oo], and vanishes only at s = U}+{a). Moreover, there is an open inter¬ 
val I"*" containing (a) such that, for any open set J C 1“'', 

lim - log UJ+ (IX G X - [ (Ts (a:) ds G ^ — inf (s). 

t^oo t t Jq } j 
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(5) The Central Limit Theorem holds. That is, for any Borel set B C 


lim io+ < a; G X 

t—fOO ' 




a-®^/2a+ 


dx 


where a+ = e" (0). 

(6) The strong law of large numbers holds. That is, for ui+-a.e. x G X, we have 

1 f* 

lim - / (Ts(x) ds = a;+(cr). 

n^oo t Jq 

(7) Let J be as in Theorem 2.4. Then e+(a) = e{a) for a G Jf^ H J. Moreover, there is an open interval 
J+ C such that /~*'(s) = I{s) for s G 1“''. 

Remark. This theorem is a refinement of Proposition 9.5 in [JPR]. We point out that parts (1) and (3) of 
that proposition are inaccurately formulated: in part (1), the interval {—5, 1 + (5) has to be replaced with 
(—5, 6), while in part (3) the interval (—((t)+ — e, {a)+ + e) has to be replaced with {{a)+ — e, (cr)+ + e). 


Finally, we have the following analogue of Theorem 2.5 on statistical properties of the dynamics under 
the limiting measure uj+. Let {f„} C K+ be an arbitrary increasing sequence going to +00 such that the 
intervals defined by (26) converge to a limiting interval j+. 

Theorem 2.7 Under the hypotheses of Theorem 2.6 the following assertions hold. 


(1) For a G , the limit 

e+{a) := lim ^et„+{a) (29) 

n—>-<50 tn 

exists and defines a real-analytic function on J^. If a does not belong to the closure of J'^, then 

lim sup — et„+{a) = 00 . 

n—foo In 


Moreover, e+{a) and tr(£'(a;)^) coincide on their common domain of definition. 

(2) The Large Deviation Principle holds in the following form. The function 

I^{s)= sup (as —e+(—a)) 

—aC:j+ 


is convex, takes values in [0, cx)] and vanishes only at s = a;+(cr). Moreover, for any open interval 
J/ C M, we have 


lim — log UJ+ 

ra—>00 tn 


X G X 



crs(x) ds G 


= - inf i+(s). 
sej 


The proof of this result is completely similar to that of Theorem 2.5, and therefore we omit it. 

Remark. Unlike in the case of the Evans-Searles symmetry, there is no a priori reason why the limiting 
intervals should be symmetric around a = |, and indeed in all cases we know where J+ can be 
computed, this property does not hold. Hence, the relation e+(a) = e+(l — a) may fail since one side may 
be finite and the other infinite, leading to the failure of the Gallavotti-Cohen symmetry /+(—s) = /+(s)+s. 
The fact that for unbounded entropy production observables the Gallavotti-Cohen symmetry may fail is 
known in the physics literature [BaCo, BGGZ, BJMS, Fa, HRS, Vil, Vi2, ZC]. In these works one can also 
find various prescriptions how the entropy production observable can be modified so that the Gallavotti- 
Cohen symmetry is restored. We shall discuss this topic in the continuation of this paper [IPS]. 
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2.6 Perturbations 

We shall consider the following type of perturbation of the reference state w. Let P be a bounded selfadjoint 
operator on JC such that D~^ + P > 0. To avoid trivialities, we assume that P is not the zero operator. Let 

+P)-i 

and let be the centered Gaussian measure with covariance . Obviously, 

DP = + Pt)-\ 


where Pt = e Pe We consider the following two cases, assuming that (G1)-(G4) hold for D. 

Case 1. P is a non-negative trace class operator such that i?P = P'd, and s - lim Pt = 0. 

>-±00 

In this case, oj^ and w are equivalent and (G1)-(G4) also hold for . Moreover, using the superscript P 
to denote the objects associated with the initial measure we easily check that 

= EP{a)=E{a), P + PC), = uj+{a), 

where we used (21) to derive the second relation. We also see that the functions e^{a) and e(a) coincide 
on J n J_^. It is possible, however, that 7 ^ J and ^ , and in fact the difference could be quite 

dramatic. Indeed, let us fix P and consider the perturbation AP for A > 0. Pick a unit vector ip such that 
Pp = ep with e > 0. 

We consider first the case of . One easily sees that for any a > 1, 

{p, + aT^P)p) < “ - A ((a - l)e - a{p, Ptp)). (30) 

m 

There exists to such that for t > to, {a — l)e — a{p,Ptp) > {a — l)e/2. Hence, for t > to and 
A > 2a/em{a — 1) the right hand side of (30) is negative which implies that a > 1 + 6t^. Thus 

= lim inf < a — 1 

t—^OO 

provided A > ^ajemioL — 1). Letting now a j, 1 we conclude that 

lim = 0 , 

A—>-oo 

and the intervals collapse to [0,1] in the limit A —> 00 . 

To deal with the case of , we set = Q^^p for a > 0 and for a < 0. A simple analysis 

yields 

(V'a.t, ((P+^)"^ - “ A|a|(e - {p, Ptp))- 

Repeating the previous argument, one shows that the length of the interval goes to zero as A — 00 , 

so that the intervals collapse to { 0 }. 


Case 2. P > 0, i?P = Pd, and Pt = P for all f G K. 

Hypotheses (G1)-(G4) again hold for D^, and we have 

P^ = (p;i + P)-i, = aP = a. 

Replacing P with AP, it is easy to see that 6^^, defined by (10), satisfies limA->.oo = 00 . Since 
{—6^^, 1 + C and C we see that the intervals and extend to 

the whole real line in the limit A —?> 00 . 
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3 Examples 


3.1 Toy model 

Suppose that the generator £ satishes £* = —£, and let G Aii be a unit vector such that the spectral 
measure for £ and tf is purely absolutely continuous. Let 


D = / + AP, 


V’ 


where Py, = (ip, ■ )(p and A > —1. Then Dt = I + where (pt = is a continuous curve of unit 

vectors converging weakly to zero as f —^ +oo. Let A± = |(|A| ± A) denote the positive/negative part 
of A. One easily verifies that (G1)-(G3) hold with m = 1 — A_, M = 1 + A+ and D± = I, so that 


6 = 


1 1 
2 +A 


Without loss of generality we may assume that (G4) holds. ^ Since (/ + XP^) ^ = I — p^P,/, for any unit 
vector ij} and any A ^ —1, we have 


D-^ PaTt =1- 


A 




1 + A 
A 


((1 - a)Py, + aPy,^), 

' ~ Ppt ) ■ 


P> Vt) 

Using the simple fact that for any two linearly independent unit vectors p, ip and all a,b G M., 


sp{aPy, + bP^) = {0} U ± 


-h 


2b{ip, p)'^ 


one easily shows that 


Sf = ^ - 


1 + A 


1 


J+ = la G 


a < 


1 + A 


4 ' x^l-{p,pt^) 2 ’ ^]xWl-[p,PtY\' 

Recalling that [p, pt) —>■ 0 as f —>■ +oo we see that for all A > —1, 6 = 6 = 5 and = (—<5+, i5+) where 

1 + A f S for A G (-1,0], 


(5+ = 


|A| 


1 + 5 for A G [0, oo). 


Furthermore, evaluating Relations (46) and (86) established below, we obtain 

et' ' 


(a) = -i log ^1 + - a) (1 - (+, . 

(1 + A)»“Ai-(*>.-/>.)"))■ 


et+(a) = -ilogfl- 


It follows that 


lim pet{a) = 
t—>-oo t 


0 for |a — A| < A + 5, 
+00 for ja—A|>A+5, 


1 I 0 for jaj < 5+, 

hm -et+{a) = < 

t I +(X) for jaj > 5+. 


Finally, one easily compute the Legendre transforms of these limiting functions, 

/(s) = (A+5)|s| - As, /+(s)=5+|s|. 

While the hrst one satishes the huctuation relation, i.e., J(s) + ^s is an even function, the second one does 
not. 


^That can be always achieved by replacing >C with A^©^, £ with £©£*, with p^ip© i/j, and setting 19(1/11 ©■i/12) = i/i2©i/’i. 
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3.2 One-dimensional crystal 


We follow [JOPP] and consider the simplest example of the one-dimensional harmonic crystal. If A C Z 
is the crystal lattice, then the phase space and Hamiltonian of the harmonic crystal are 

= {{p,q) = ({p„}„gA,{gn}nGA) \ Pn,qn G R}, 


HA{p,q) 

riGA ^ 


gn , (gn-gn-l)^ 
2 ^ 2 


where we set qk = 0 for k ^ A (Dirichlet boundary conditions). The Hamilton equation of motions are 




where 


Ca = 


0 -JA 
1a 0 


jA is the restriction of the finite difference operator 


{jq)n — gn-t-1 gn—1 (31) 

to with Dirichlet boundary condition and 1 a the identity on (which we shall later identify with 
the projection > K^). Clearly, for all A, Ja is a bounded selfadjoint operator on ^^(A) satisfying 
1 < JA < 5. 

To fit this model into our abstract framework, we set Ta = Ax Z 2 , Ta = R'""' = © R^ with the 

weight sequence I = where In^i = ca(1 + and ca is a normalization constant. One 

easily verifies that Ca^a; ^ ^Ai dynamics of the harmonic crystal is described by the group . 

Let /lA be the self-adjoint operator on /Ca = ^b(^) ® ^b(^) associated to the quadratic form 2Ha- Energy 
conservation implies C\hA -f /ia>Ca = 0. Equivalently, the operator La defined by 

T — r — ( ^ ~jA^\ 

La — hj^ Ea^a — 0 j ’ 

is skew-adjoint. Since 1 < Ha < 5, this implies in particular that the group is uniformly bounded on 

/Ca. 

Our starting point is harmonic crystal on A = Z and in this case we drop the subscript A. Eor our purposes 
we will view this crystal as consisting of three parts, the left, central, and right, specified by 


Ae = (-00,-1], Ac = { 0 }, Ar = [ 1 ,00). 


In what follows we, adopt the shorthands Ha/. = Hi, Hai = hi, j'a^ = ji, etc. Clearly 


X = © Xc © Xr, JC — K-l © JCc © /Cr, 

where JCs = ^b(^s) ® ("bI^s) s = I, c, r, and 

H = Hq + Vi + 1A, 


where 

Ho = Hi+ H, + Hr 
and Vi{p,q) = -qoq-i, Vr{p,q) = -gogi- 

The reference state uj is the centered Gaussian measure with covariance 

D = Di® Dr® Dr, 
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where 




s = i,c, r, 


Is is the identity on £^(As), and Tg > 0. Thus, initially the left/right part of the crystal are in thermal 
equilibrium at temperature T^/r- The Hamiltonian Vnjr couples the left/right part of the crystal to the 
oscillator located at the site n = 0 and this allows for the transfer of the energy/entropy between these two 
parts. The entropic fluctuation theorems for this particular Gaussian dynamical system concern statistics 
of the energy/entropy flow between the left and right parts of the crystal. 


Hypothesis (G1)-(G4) are easily verified following the arguments of Chapter 1 in the lecture notes [JOPP] 
and one finds that 


uj+(a) 


XTt-Trf 
' T,Tr 


where k = (v^ — l)/27r, and 


e{a) 


—ttlog 



[Tj - Trf 
T^Tr 


a{l 



(32) 


Note that e{a) is finite on the interval Jo = (—(5o, 1 + ^o), where 


_ mm{Ti,Tr) 
~ \Ti-Tr\ 


(33) 


and takes the value +oo outside the interval Jq. Note also that 6o can take any value in (0, oo) for appro¬ 
priate choices of Tg,Tr G (0, oo). The measure i/ in Part (2) of Theorem 2.4 is 


1/ = kD-Ss + 


where 'Da is the Dirac measure centered at a. 

We finish this section with several remarks. 

Remark 1. The intervals J, can be strictly smaller then Jq. To see this, fix Tc, 6o, a > 1, and set 
Tr = (1 + to ensure Relation (33). Let ip G JChe such that {ip, hcp) = 1. One has 

{p, + aTt)p) = X! ^ ’ 

where pt = e~^^p. Since the skew-adjoint operator L has purely absolutely continuous spectrum and 
is compact, there exists to > 0 such that 

{pt,hcPt) = 

Za 

for all t > to- Moreover, since the Hamiltonian flow is uniformly bounded there exists a constant C such 
that 

1 Ol 

((1 - a){p,hi/rp) + a{pt,heirPt)) <C—. 

J-e/r 

Summing up, if Tg > ACToaj{a — 1), then 

{p, [D-^ + aTt)p) < ^ + 20^ < 0, 

for all t > to and hence S < a. Thus, in the limit Tg —g oo the interval J collapses to [0,1]. In a similar 
way one can show that in the same limit the interval collapses to {0}. On the other hand, arguing as in 
the Case 2 of Section 2.6, one can always take Tg/r,Tc -G 0 in such a way that in this limit the intervals J, 
J'^ extend to the whole real line. 
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Remark 2. Somewhat surprisingly, even in the simplest example of the harmonic crystal discussed in 
this section, it appears difficult to effectively estimate the location of the intervals J, outside of the 
perturbative regimes. In particular, the subtleties regarding the location of these sets were overlooked in 
Sections 1.11, 1.14 and 1.15 of the lecture notes [.I OPP] . These difficulties raise many interesting questions 
and we leave the complete analysis of these aspects as an open problem. 

Remark 3. An interesting question is whether one can find P such that for the perturbed reference state 
as defined in Section 2.6 one has J = Jo. That can be done as follows. Set /3s = l/T/j, suppose that 
Pr > Pi and let 

^ ({^r - /3c)lc 0 

0 (/3r + 2j3(, - 3/3c)jc + fitvt + PrVr 

where vi/j. denotes the selfadjoint operator associated with the quadratic form One easily checks 

that 

= {prh- Xh‘f '>)-\ 


where X = Pr — Pi > Q, 


,{N) _ /lAfUA, 

“(v 0 



and denotes the restriction of the operator (31) to with Neumann boundary condition. We are 

concerned with the interval 


Jf = {a G K I {D^)-^ + aT[ > 0}. 


Since 


(of )-^ = Prh - Xe-^^" ^1/2^ 


a simple computation gives 

(O^)-i + aTf = - (1 - 

and hence 


jf = {a G K I/3r/2f > (1 — a)/i + 


Since PrjX = 1 + (5o and 

0 < hf'' < h, 

we have that for all t, 

{-5o, 1 + <5o) C jf. 

Thus, limt_>oo and = Jo- 

Remark 4. In contrast to Remark 3, we do not know whether there exists P such that for the perturbed 
reference state ui^ one has = Jq. 

Remark 5. In the equilibrium case Ti = Tr = T have 02 +(cr) = 0, and one may naively expect 
that a does not fluctuate with respect to uj and w+, i.e., that e(a) = e+(a) = 0 for all a, and that 
I{s) = I+{s) = 00 if s 7 ^ 0. If one also takes = T and the perturbed reference state described in 
Remark 3, then cr = 0, and the above expectation is obviously correct. On the other hand, for the reference 
state determined by D, in the high-temperature regime T —> 00, fixed, the interval J collapses to [0,1] 
while the interval collapses to {0}. Hence, in this regime, the rate functions I{s) and I'^{s) are linear 
for s < 0 and s > 0 , with the slopes of the linear parts determined by the end points of the finite intervals 
J and J+, and the entropy production observable has non-trivial fluctuations. 

Remark 6. The scattering theory arguments of [JOPP] that lead to the derivation of the formula (32) extend 
to the case of inhomogeneous one-dimensional harmonic crystal with Hamiltonian 


H^ip,q) = J2(Y 

neA ^ 


2 


i(^n 


- g»-l)^ 
2 ^ 
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where w„ and are positive numbers satisfying 

C~^ < Kn < C for all n G Z, 
and C > 1 is a constant. In this case the operator j is the Jacobi matrix 

(j9)n — (^n tvyi -f ^nQn—1 ^n+l^n+1; ^ G Z. 

One easily verifies that Hypotheses (Gl), (G2), and (G4) hold. If j has absolutely continuous spectrum 
(considered as a self-adjoint operator on ("^(Z)), then (G3) also holds. Moreover, a;+(cr) and e(a) can be 
computed in closed form in terms of the scattering data of the pair (j, jo), where jo = jr© jc©jV (for related 
computations in the context of open quasi-free quantum systems we refer the reader to [JLP, JOPP, Lan]). 
The formulas for uj+{a) and e(a) involve the scattering matrix of the pair (j, jo)^ estimating the 
location of the intervals J, J’*’ is difficult. However, the interesting aspect of the formula for e(a) is that it 
allows to express the measure i/ in Part (2) of Theorem 2.4 in terms of the scattering data. The mathematical 
and physical significance of this representation remain to be studied in the future. Finally, the scattering 
methods can be extended to treat an arbitrary number of infinite harmonic reservoirs coupled to a finite 
harmonic system. The discussion of such extensions is beyond the scope of this paper. 


4 Proofs 

4.1 An auxiliary lemma 

Using the notation and conventions of Section 2.1, we have the following simple result. 

Lemma 4.1 (1) If A = A* G 'T, then the quadratic form <Ia(a) = (x,Ax) has a unique 

extension to an element of (X, doJo) with a norm satisfying ||fj'. 4 ||i < ||U|| ||^||i- Moreover, 

[ qA{x)duJD{x) = tr{DA). (34) 


(2) Let K 9 f i-G At = G T be differentiable at t = to and let Atg be its derivative. Then the map 
ffi 9 f I—>■ qAt G L^(X, dwo) is differentiable at t = to and 


dt 


= 9 ©, 


t=to 


(3) If 1 does not belong to the spectrum of A, then the function T 9 i-9 F{X) = det(J — X) is 

differentiable at X = A and its derivative is given by 


(D^F)(X) = -F{A) tr((/ - A)-^X). 


(35) 


Proof. 

Part (1) By Eq. (4), the function x i-9 ^y{x) = {y,x) belongs to L^(X, dw£)) for y G X;*. Moreover, 
Fubini’s theorem yields the estimate 

W^yWl = 51 [ XiXjdujDix) = D.^ytyj = {y,Dy) < \\D\\ ||j/f, (36) 

i,jer i,jer 

which implies that the linear map y i-G has a unique extension $ ; ^R(r) —>■ Lf{X, dujo), such that 
||$|| < ||i9||V2. 

^In the case of harmonic crystal considered in this section, j is a discrete Laplacian and the absolute values of the entries of the 
scattering matrix of the pair (j, jo) either O’s or I’s. For this reason the formula ( 32 ) for e(Q:) has a particularly simple form. 
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A self-adjoint A G T has a spectral representation A = '^k‘Pki‘Pk, ■), where the Uk are the eigenvalues 
of A and the corresponding eigenvectors ipk form an orthonormal basis of ^^(r). It follows that qA{x) = 
Thk from which we conclude that qa extends to an element of dujo) with 

\\qA\\l<Y.M\\<i>^Jl<J2\<^k\\\D\\ = \\D\\ Pill. 

k k 

The last equality in Eq. (36) yields 

/ qA{x)du}D{x) = PV 3 JI 2 = '^ak{(pk,D<fk) = tr{AD), 

k k 

which proves Identity (34). 

Part (2) It follows from Part (1) that the linear map T B A i-A qA G L^{X, dwn) is bounded and hence C^. 
Part (3) Using a well known property of the determinant (see Theorem 3.5 in [Si]), we can write 

F{A + X) = det(/ - (A + X)) = det((/ - A){I -{I- A)-p) 

= det(/ - A) det(/ -{I- A)-P) 

= F{A) det(/- (/-A)-P). 

To evaluate the second factor on the right-hand side of this identity, we apply the formula 

OO 

det(/ + Q) = l + y]tr(Q^'=), 

fc=i 

where denotes the /c-th antisymmetric tensor power of Q (see [Si]). Since ||(5^^||i < (A:!)“^||Q||i, 
one has the estimate 

I det(/ + Q) - 1 - tr(g)| < elKSIIi _ ^ _ ||g||^ < _|jQ||2. 

It follows that 

det(J - (7 - A)-P) = 1 - tr((7 - A)-p) + 0(\\X\\\), 
as X 0 in T. Thus, we can conclude that 

F{A + X) - F{A) = -F{A) tr((J - A)-p) + 0(|1X||2), 
and the result follows. □ 


4.2 Proof of Proposition 2.1 


Part (1) Up to the constant tr(£)<;) (which is well defined since G F), a is given by the quadratic form q^ 
which is in L^iX, dw) by Lemma 4.1 (1). For x G Xi, i.e., w-a.e. x G X, one has 


crt(x) - as{x) = \{x,{ 




whence, setting q = and applying again Lemma 4.1 (1), it follows that 


Ikt - (J,\\Li(X,dF < ^\\D\\ Ikt - Cslll. 

Thus, it suffices to show that the function f i—> G T is continuous. This immediately follows from 
the norm continuity of the group the fact that G T, and the well-known trace inequality pB||i < 
Pll Pill. We note, in particular, that 

lk*|Ui(ie.d.)<PII(l + ||kkP)lklli forfGM. 
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Part (2) From Eq. (5), we deduce that 

Ctlo; = ^ logdet(/ + DTt) (37) 

Now note that Tt = — D~^ satisfies the cocycle relation 

Tt+s=Tt+e-*^'Tse-*^. (38) 

It thus follows from Assumption (Gl) that the function t Tt G T is everywhere differentiable and that 
its derivative is given by 

ft = (39) 

Lemma 4.1 (3) and the chain rule imply that the first term on the right-hand side of (37) is a differentiable 
function of t. Using Eq. (35), an elementary calculation shows that 


^||^logdet(/ + DTt) 


-tr(Z70- 


Applying Lemma 4.1 (2) to the second term on the right-hand side of Eq. (37), one further gets 


Summing up, we have shown that 


' 2 df 
d „ 


<lTt = =q,o, 


= cr-t, tGR. 

Since the function t i—>■ a-t G dcj) is continuous by Lemma 4.1 (1), and 

Riemann’s integral to write 

Lt\cj= / cr_sds. 

Jo 

The fact that, for w-almost every x G X, one has 

^ujt\uj{x)= / fj_s{x)ds, 

^0 


0 , we can use 

(40) 

(41) 


follows from Theorem 3.4.2 in [HP]. 

Part (3) Prom the cocycle relation 

~~i~ °‘P *) ( 42 ) 


we infer 


and hence 



a-t e 



— 1 — sa) o (j)' 


t 

dw ’ 


ICsl dw 



1 — sa 


dw < 




^UJs \<JJ 



^U)s\UJ 


dw. 


To prove that Relation (9) holds in (T, dw), it suffices to show that both terms on the right-hand side of 
this inequality vanish in the limit s —0. 

To estimate the first term we note that the inequality - 1 - ^ > 0 (which holds for £ £ K) combined 
with Eq. (34) and (37) implies 


J - 1 dw = 1^ - 1 - y £<^,|^dw^ 


(tr(Z9r,)-logdet(/ + Z7r,)) 
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By Assumption (Gl), the map s i—Tg is differentiable in T at s = 0. Since Tq = 0, we can write 


lim — I |e “""I" - 


du=- 


— (tr(OTg) -logdet(/ + OTg)) 
ds 


s =0 


s-J-O |S 

Using Lemma 4.1 (3) and the chain rule, we get 

^ (tr(DTg) - logdet(/ + DTs)) = tr(OTo) - tr(OTo) = 0. 
ds s=o 

To deal with the second term, we use Eq. (40), Eubini’s theorem and Lemma 4.1 (1) to write 


pi p 

< / 


Jo 

Jx 


) du 




Ik-™ - ‘^llidu, 


and since the map s i —q is continuous in T, the dominated convergence theorem yields 


lim 

s—>-0 


k-™ - ?|iidu = 0. 


Part (4) Relation (7) implies that 


wkcr) = w(CTt) = / do; - ti{D‘;), 
Jx 


and formula (34) yields 


ujt{a) = - <;)) = tr(<j(A - D)). 

Part (5) Starting from Definition (2) and using the cocycle relation (42), we obtain 

Ent(wt|w) = - / iujt\ujdu}t= / C_^|^dw. 

Jx Jx 


Eq. (41) and Fubini’s theorem further yield 


Ent(a;i|a;) = / / ct-s dsdcj = — / / (7^ dsdu; = — / a;s((7 )ds. 

Jx Jo Jx Jo Jo 


4.3 Proof of Proposition 2.2 

Part (1) We have to show that uj+, the Gaussian measure of covariance D^, is the weak limit of the 
net {wiloQ. Since the cylinders form a convergence determining class for Borel measures on X (see 
Example 2.4 in [Bill]), it suffices to show that = uj+{Ci{B)) holds for any finite 

subset / C r and any Borel set i? C K^. By Hypotheses (G2)-(G3), one has limt_>.oo Dtj = 7?+,/ and 




< e 2 “ , 


for all X € It follows that limt_).oo as well as limt_>oo det(27ri7tj) = det(27ri7+j) so 

that 


lim 


1 


r / e dx = ^ / e da;, 

Jb A/det(27rD+j) Jb 


y/det{2TTDtj) Jb 

holds by the dominated convergence theorem. The same argument applies to uj_. 
Part (2) Eollows directly from Lemma 4.1 (1) and Proposition 2.1 (4). 
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4.4 Proof of Proposition 2.3 

Part (1) Let us note that a G Jt if and only if 

D-^ + - D-^) > 0. (43) 

It follows that Jt is open. For 0 G [0,1], we can write 

D-^ + eaie -*^'= 9 (d"! + - D"!)) + (1 - 9)D~\ 

whence a G Jt ^ 9a G Jt and we can conclude that Jt is an interval. Multiplying (43) by -d from the left 
and the right and using the relations = r?* = we obtain 

+ >0, (44) 

whence we see that a G J-t- By symmetry, we conclude that J_i = Jt- Furthermore, multiplying (44) 
by and from the left and the right, respectively, we obtain 

aD-^ + (1 - > 0 . 

It follows that 1 — a G Jt, and by symmetry, we conclude that a G Jt if and only if 1 — a G Jt- Thus, Jt 
is an open interval symmetric around a = ^. 

Part (2) For any bounded operator C > 0 on ^^(r) and for any a,t G M such that C~^ + aTt > 0, 
formulas (5) and (37) allow us to write 


dwc 


(det(I + DTt))" 
det(I + aCTt) 


(45) 


By definition D ^ + aTt > 0 for a G {—St, 1 + St). Taking C = 77 in (45) and integrating over X, one 
easily checks that 

ei(a) = ^ logdet(/ + DTt) - ^ logdet(/ + aDTt) (46) 

for all f G K and a G (—St, 1 + St). The first term on the right-hand side of this identity is linear in a 
and hence entire analytic."^ The determinant in the second term is also an entire function of a, and its 
logarithm is analytic on the set where the operator I + aDTt is invertible; see Section IV. 1 in [GK]. 
Writing I + aDTt = D{D~^ + aTt), we see that I + aDTt is invertible for a G Jt- Furthermore, since 

I + aDTt = a77i/2(a-i/ -p D^/^TtD^/^)D-^/^, 

and the operator D^/'^TtD^^'^ is self-adjoint, we conclude that I -f aDTt is invertible for a G C \ K. 
Hence, the function et{a) is analytic in the cut plane C+ UC_ U Jt- Its convexity is a well-known property 
of Renyi’s relative entropy and follows from Holder’s inequality applied to Eq. (12), and relations (13) are 
easy to check by a direct computation. 

It remains to prove that et{a) = -foo for a ^ Jj. To this end, we first note that the spectrum of D~^ 
is contained in the interval [M~^,m~^] and that the operator aTt is compact. By the Weyl theorem on 
essential spectrum, it follows that the intersection of the spectrum of the self-adjoint operator D~^ -f aTt 
with the complement of [M~^,m~^] consists of isolated eigenvalues. Thus, if a ^ Jt, then there are 
hnitely many orthonormal vectors numbers \j > 0, and an operator B > cl with c > 0 such that 


D ^ + aTt — — 'y ^ Xj{tpj, -f B- 

7=1 


^We shall see in the proof of Theorem 2.4 that it is in fact identically equal to zero. 
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It follows that 

= {det{I + DTt)^^^ J expj ^ ^ (47) 

Since B — D~^ G T and D~^ + B > 0, we conclude from (5) that coincides, up 

to a numerical factor C > 0, with a centered Gaussian measure whose covariance operator is equal to 
£>' := (L>~^ + B)~^. Hence, we can rewrite (47) in the form 


f exp|i^Aj|(v3j,a:)|2|u;_D/(dx). 


Since the support of ujui coincides with the entire space, this integral is infinite. 

Part (3) Using the cocycle relation (42), we can write ^ 

et(l — a) — log = loga;t(e““^‘"‘i‘^) 

= log) = logtLi(e“^“-*i“) = e_t(a). 

Now note that, by (G4), the measure w is invariant under 'd, whence we conclude that = oJt o'd and 
^wt\uj o t? = It follows that e_t(a) = et{a). Combining this with the above relation, we obtain the 

Evans-Searles symmetry. 


4.5 Proof of Theorem 2.4 

Part (1) We first prove the existence of limit (14). Let us set 

Dt{a) = {{I-a)D-^+aD-^)-^ (48) 

and recall that et{a) can be written in the form (46). Using Relations (35), (39), Lemma 4.1 (3) and the 
chain rule we obtain 

^ logdet(/ + aL>Tt) = tv{{I+ aDTt)~^aDft) =-2a tr(A(a)<t_t) =-2a tr(L)_t(l-a)^). (49) 

In particular, for a = 1 the derivative is equal to zero for any f S M, whence we conclude that the first term 
in (46) is identically equal to zero. Let us now fix a S J and choose fo > 0 so large that a G Jt for t > to. 
It follows from (46) and (49) that 

1 1 2 q( 

-et{a) =-etoia) - — tr(i9_s(l - a)^)ds. (50) 

t t t 7to 


By Assumption (G3) 

s-limi9_s(l — a) = D_{1 — a) := {aD~^ + (1 ~ ol)DZ^) 
and since 1 ; is trace class, it follows that 

lim tr(i9s(l — a)^) = tr(i9_(l — a)<t). 

Combining this with (50), we conclude that for a G J, 

lim -et(a) = —2atr(iA_(l — a)<j). (51) 

t-f+oo t ^ ' 

^Note that this computation does not use (G 4 ). 
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Once the existence of limit is known, we can easily obtain the required properties of e{a). The convex¬ 
ity of e{a) and the first and last relations in (15) follow immediately from the corresponding properties 
of et{a). Furthermore, it follows from (40) and the invariance of u) under d that 

e't(O) = f ^ujt\oj{x)u}{dx) = 

In view of Part (2), the limit e{a) is analytic on its domain of definition. By Theorem 25.7 in [Rock], 


/ / (T_s(a:) ds w(da;) = — / / CTs(x) ds a;(dx). 

Jx Jo Jx Jo 


lim - 

>-oo t 


e'i(a) = e'(a), 


for a G J_. Using Fubini’s theorem and Part (2) of Proposition 2.2, we derive 


1 If 

e'(0) = lim -6^(0) = — lim - / a;(CTs)ds = —tu+(cr) = —tr(<ji7+). 
t —>-00 t t —>-00 t Jq 

The third relation in (15) now follows from the fourth one. 

Part (2) The analyticity of e(a) follows from Relation (51). We now prove (16). 

1/2 1/2 

Let fj, be the spectral measure of Q for the linear functional induced by the trace class operator D_ <;D_ . 
In other words, /r is the signed Borel measure such that 


J = tr(/(Q)£)y\£)y^), 


(52) 


for any bounded continuous function / : M —C. By Eq. (20), the measure /r has its support in the interval 
[—(5 ,(l + d)“^]. One easily checks that 


defines a continuous linear functional on the Free bet space C'o(ffi) of compactly supported continuous 
functions / : M —>^ C. By the Riesz representation theorem (see Chapter 2 in [Rud]), it follows that there 
exists a signed Borel measure ly, with support on (—oo, —<5] U [1 + d, c»), such that 

I /(r)i.(dr) = J /((?-i)g-V(dg) (53) 


A standard argument based on the monotone class technique shows that (53) remains valid for any bounded 
measurable function /. Decomposing the measures p and i/ into their positive and negative parts, we easily 
deduce from (53) that 

J /(r)|H(dr)= J f{q-^)\qn^\{dq), 

for all bounded continuous /. In particular, taking /(r) = - outside a small neighborhood of zero and 
using (52), we derive 

J = j < 00. 

Recalling relation (22) (which will be established below) and using (53) with /(r) = — log(l — ar~^) on 
the support of i/, we obtain 

e{a) = = — J ag{aq)ii{dq) 

= — J q~^ log(l — aq)g{dq) = — J log(l — ar~^)i'{dr). 

This relation coincides with (16). 
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To prove the uniqueness, let i^i, V 2 be two signed Borel measures with support in K \ J, satisfying 
/ < oo, fc = 1, 2, and such that 


J log(l — ar ^)i^i(^dr) = J log(l — ar ^)i/ 2 (dr) 


for a G J^. Differentiating, we derive that 


/ 


dt^i(r) 
r — a 



dr^W 

r — a 


(54) 


for a S J. By analytic continuation (54) holds for all a G C+ U C_. Since the linear span of the 
set of functions {(r — a)~^ \ a G C+ U C_} is dense in (7o(K), (54) yields that for any / G (7o(K), 
/ fdvi — J /dz/ 2 - Hence vi = V 2 . 

Part (3) The fact that / is a convex function taking values in [0, + 00 ] follows immediately from the def¬ 
inition. The relation e'(0) = oj-{a) = —a;+(cr) and the regularity of e imply that I vanishes only at 
s = uj+{a). The validity of (17) is a straightforward consequence of the last relation in (15). Let us 
prove (18). 

Consider the following family of random variables {St}jg[o defined on the probability space (X, T, uj) 


St = 



(Ts ds. 


By Proposition 2.1 (2) and the symmetry relations to = lu o-d and a o d = —a, we have 


et(a) =loga;(e“'^-*i“) = logw(e“-/'o= logu;(e-“-/'o= loga;(e-“*^*), 


so that et{—a) is the cumulant generating function of the family {St}tg[o,oo)- Applying a local version of 
the Gartner-Ellis theorem (see Theorem 4.65 in [JOPP]), we conclude that (18) holds with 


£ = min(— a;+(tT) — d'^e{—d), —uj+{a) + d e(l + S)) = min(e'(0) — d'^e{—S), d e(l -f 5) — e'(l)) 

where d^e{a) denotes the right/left derivative of e{a). The fact that £ > 0 follows from the convexity and 
analyticity of e(a). 

Part (4) As was shown above, et{—a) is the cumulant generating function of {Ej}. Therefore, by Bryc’s 
lemma (see [Br] or Section 4.8.4 in [JOPP]), the CLT will be established if we prove that et{a) extends 
analytically to a disc Vg = {a G £.\\a\ <£} and satisfies the estimate 


sup -\et{a)\<oo, 

t'>to,CtG'De ^ 


(55) 


for some > 0. The analyticity was established in Part (2) of Proposition 2.3. Using the representa¬ 
tion (50), one easily sees that in order to prove (55) it suffices to show that 


sup ||A(a)||<oo. 
11 —ck| 


(56) 


An elementary analysis shows that Assumption (G2) implies the lower bound 

(1 _ + „D-. > 1 (a + 1 _ l„ _ 11), 


M M + m 


for t,s gM. and a G [—5,1 -f <)]. Since for z £ 


Re {{1- z)D-^ + zD^^) = {l-Rez)D-^+RezD:;^\ 


we have the upper bound 


a, , 1 M M + m , ^ ^ , ,,, 

((l-z)D^ + zDf ) 2I) 


-1 


(57) 


(58) 
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for s,t gM. and z in the strip {z S C | Re z G (—15,1 + (5)}. Thus, the required estimate (56) holds provided 
e < 6. 

Part (5) We first note that the differentiability of e{a) at zero and a local version of Theorems II.6.3 in [El] 
(which holds with identical proof) implies that, for any £ > 0 and any integer n > 1, 

uj{{x GX \ |E„ - a;+(tT)| > £}) < 


where a{£) > 0 does not depend on n. By Theorems II.6.4 in [El], it follows that 

1 r 

lim — / CTs(a;) ds = a;+((T) 
n-foo n Jq 

for w-a.e. x G X. Suppose now we have shown the following inequality for some r < 1 

/*n+i 


sup / (Ts(a:)ds 
0<t<l Jn 


< (n + l)" 


for n > no(x), 


(59) 


(60) 


where no{x) > 0 is an integer that is finite for uj-a.e. x G X. In this case, we can write 


If*,, 1 ^ 

1 


1 

r . ^ 

- / as{x)ds -/ crs(x)ds 

< - 

/ crs(x)ds 

+ ^ 

/ crs(x)ds 

t Jo n Jo 

n 

Jn 


Jo 


where n is the integer part of t and t = t — n. It follows from (60) that the hrst term on the right-hand 
side goes to zero for a.e. x G X, and the second goes to zero in view of (59). Combining this with (59), we 
obtain (19). Thus, it remains to establish (60). 

Let us fix an arbitrary r G (0,1) and denote by ^„(x) the expression on the left-hand side of (60). In view 
of the hrst relation in (11), we have 


^„(x) = sup 

0<i<l 


pn-\-t 

/ ^e®'^x)ds 

J n 



rn+t 

= sup (x,>t„,tx) , 

1 ds 

0<i<l 

Jn 


Suppose we have constructed a sequence {Bn} of selfadjoint elements of T such that, for any n > 0, 

sup |(x,<J„,tx)| < {x,BnX), \\Bn\\l < C, (61) 

0<t<l 


where C > 0 does not depend on n. In this case, introducing the events An = {x G X | ^n{x) > (n-f 1)’’}, 
for sufficiently small £ > 0, we can write 

uj[An) < (det(J - 2£DR„))■^^^ (62) 

where we used the fact that the Gaussian measures on X with covariance operators — 2eBn)~^ 

and D are equivalent, with the corresponding density given by (see (5)) 

Ad'|d(^) = (det(/- 


In view of the second inequality in (61), the determinant in (62) is bounded from below by a positive 
number not depending on n > 0 for sufficiently small £ > 0. Thus, the series converges, and 

by the Borel-Cantelli lemma, inequality (60) holds with an almost surely hnite integer no(x). 

We now prove (61). Erom Assumption (G2) we derive 


so that the uniform bound 


lle‘^|| < 


,t£* 

> me‘ 

'M 

\ 1/2 

m 

1 ’ 


(63) 
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holds. Since G T is selfadjoint, one has |(x, ^a:)| < (x, |^|a;) for all x G 1C. Hence 


where 


sup < / |(e®'^x, <je'’'^a:)|ds < / |<;|e®^a;)ds = (a;, i3„x), 

J n J n 

rn-\-l 

J n 


Br, = 


is a self-adjoint element of T such that 


The proof of Theorem 2.4 is complete. 


M 

m 


4.6 Proof of Theorem 2.5 


Part (1) Let {s„} be an arbitrary sequence converging to 6. Recall that D ^ > 0 for a G Jg^. 

Multiplying this inequality by from the right and by e®"^ from the left, we obtain 

(1 - > 0, 

for any a G Js„. Invoking Assumptions (G2)-(G3), we can pass to the limit in the last inequality to get 

(1 - a)DZ^ +aD-^ > 0, 

for any a G J. Taking a = 1 + S and a = —6 and performing some simple estimation, we obtain 
inequality (20). Furthermore, it follows from (20) that aQ < 1 for a G (—S, 1 + 6), whence we conclude 
that the operator function (21) is analytic in the cut plane C+ U C_ U {—5 ,1 + J). 

Part (2) We first prove the existence of the limit in (22). To this end, we shall apply Vitali’s convergence 
theorem to the sequence of functions 

hnia) = ^et^{a), n > 1, a G Jt„. 

By the very definition of 6, for any e > 0 there is such that, for all n > N^, the function is analytic 
in the cut plane C_ U C+ U Jg where 

Je — (— ^ + £, 1 + (5 — e) C Jt„ ■ 


By the proof of Part (4) of Theorem 2.4 (more precisely Eq. (58)), the functions /i„ are uniformly bounded 
in any disk or radius less than 6 around a = 0. By the Cauchy estimate, the same is true of their derivatives 
h' 

"'n- 

Let iTo be the compact subset of (C_UC+U Je)\{ 0} described on the left of Figure 1. From Definition (48) 
we infer 

A„(a) = + aQn)-^D^^^ = zD^/^{z - z = 

a 

where is a selfadjoint element of T. By definition, a G Jt„ iff I + aQn > 0, i.e., 

sp(Qn) C (-(1 + C (-(l + 5-e)"\((5-£:)"^) (64) 


for all n > N^. Since the function a i-G- z = —1/a maps Kq to a set which is uniformly separated from 
sp(Q„) (see Figure 1), it follows from the spectral theorem that 


sup ||A„(a)|j < ||L>|| 

n>N^ 

aGKo 


n>N, dlSt(2:,Sp((y„)) 
-z-iGifo 
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Im a Iin z 




Figure 1: A compact region Kq C (C_ U C+ U J^) \ {0} and its image under the map a i—>■ z = —1/a. 
The thick lines in the a-plane are the cuts M \ Jg. By Eq. (64), if n > Ng, then the spectrum of Qn lies 
inside the thick line of the z-plane. 


Applying Lemma 4.1 (3) to Eq. 46 (recall that the hrst term on the right hand side of the latter vanishes) 
and integrating Eq. (39) to express Tt^ we obtain 

1 

h/{a) = - —tr(A„(a)Tt„) = / tr(£)t^(a)<r_st„)ds. 

Jo 

The bound (63) further yields 

|/^;(a)|<-|k||i||A.(a)||, 

m 

and the previous estimate allows us to conclude that the sequence is uniformly bounded in Kq. 

Summing up, we have shown that {h'„}n>Ne is uniformly bounded on any compact subset of C_ UC+U Jg 
and since /i„(0) = 0, the same is true of the sequence {hn}n>Ne- By Part (1) of Theorem 2.4, the 
sequence {/i„(a)} converges for a G J. By Vitali’s theorem (see Section LA.12 in [GR]), we conclude 
that the sequence {/i„} converges uniformly on any compact subset of C_ U C_|_ U Jg, and the limit is an 
analytic function on it. Since e > 0 was arbitrary, we see that the middle term in (22) is well defined for 
any a G C_ U C+ U J and is an analytic function on this domain. 

To prove the second equality in (22), it suffices to establish it for a G J, because both left- and right-hand 
sides are analytic functions on C_ U C+ U J. The lower bound (57) shows that Dt{a) is bounded and 
strictly positive for all f G K and a G (—15,1 -f i5). It follows from Eq. (37) and Lemma 4.1 (1) that 
dujDt{a))- Moreover, Eq. (45) shows that for / G L^{X, dujDt{a)) 






(65) 


Using this relation with / = integrating the identity 

Jo 

against w, and applying Fubini’s theorem, we obtain 


= 1 + / d7. 

Jo 


Rez 
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Resolving this integral equation (which reduces to a linear differential equation) for a i—>■ we 

derive 

Taking the logarithm, dividing by t, and using (41), we obtain 

-et{a) = - W£,^( 7 )(Ct|<,;)d 7 = T / / ^£,^( 7 )(^-s) dsd 7 = / / ujD,(j){cr-ts)dsd'y. (66) 

^ ^ Jo ^ Jo Jo Jo Jo 

It follows from (34) and the first relation in (11) that 

= tr(A(7)<^-t^) = tr(e-‘®^A(7)e”‘"^\) = tr (((1 - j)Dzls + ■ 

Combining this with Hypothesis (G3) and a continuity property of the trace, we derive 

= tr(D^<j) = w;o^(cr) for 7 e (-<5,1 + 5), s G (0,1), 



where we set = ((1 — 'y)D_^ + 'yD_^^) The bound (58) allows us to apply the dominated conver¬ 
gence theorem to Eq. ( 66 ), and conclude that 

^ pCt /'I nOi 

e{a) = lim -et(a) = / / (tr) dsd 7 = / tr d 7 , a G (—5,1-I-<5). (67) 

t Jo Jo '' Jo 

Writing — 'yQ)~^D]/^, we further get 

e{a) = f 
Jo 

and performing the integral yields Eq. (22) for a G (—<5,1-1-5). 

Einally, to prove (23), it suffices to note that if a does not belong to the closure of J then, for infinitely 
many n > 1, a ^ Jt„ and by Proposition 2.3 (2), et„ (a) = -foo. 

Part (3) The required properties of the rate function I follow from (15) and elementary properties of the 
Legendre transform. Thus, we shall only prove (25). In doing so, we shall assume that the interval J is 
finite; in the opposite case, the result follows immediately from the Gartner-Ellis theorem; see Section 4.5.3 
in [DZ]. Moreover, we shall consider only the non-degenerate situation in which oj+(<j) > 0. The analysis 
of the case uj+ (tr) = 0 is similar and easier. 

Let us extend e{a) to the endpoints of the interval J by the relation 

e{a) = limsup -et(a), a G {—5,1 -f 5}. 

t —>+00 t 


Since the extended function e is convex and, hence, continuous at any point where it is finite, the Legendre 
transform of e{—a) coincides with I defined by (24). In view of a well-known result on the large deviation 
upper bound (e.g., see Theorem 4.5.3 in [DZ]), the following inequality holds for any closed subset F C'R: 


lim sup 

n—>-00 


logo; 


a; G X 



as{x) ds G F 


< — inf i(s). 

sGF 


Since I is also continuous, this upper bound easily implies that (23) will be established if we prove the 
inequality 

liminf — log w ( I a; G X — / crs{a;) ds G OI ) > — inf/(s), (68) 

n^oo tn V L tn Jo ) J 

where O C K is an arbitrary open set. A standard argument shows that it suffices to prove ( 68 ) for any 
open interval J' C M. Let us set 

s~ = — lim_e'(a), = — lim_e'(a). 
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In view of the local version of the Gartner-Ellis theorem (see Theorem 4.65 in^ [JOPP]), relation (25) is 
true for any interval J C (s“, s’*"). Thus, it suffices to consider the case when J = = (s — e, s + e), 

where ±(s — s±) > 0. The proof of ( 68 ) is divided into several steps. 

Step 1: Reduction. We first show that the required inequality will be established if we prove that, for any 
s S K satisfying the inequality ±(.§ — s±) > 0 and any e > 0 , 

liminf — loge)) >—/(s±e), (69) 

ra —>00 ' ' 

where i3„(s, e) = {a; S X | + s| < e}. Indeed, we have 

{ —(l + 5)s —e“ fors<s“, 

. \ ^ (70) 

(5s — e+ for s > s+, 

where e~ (respectively, e’*') is the limit of e(a) as a f 1 + (respectively, a i —6). In particular, the rate 
function I is everywhere hnite and continuous. It follows from (69) and inequality ( 68 ) with J' C (s“, s+) 
that 


lim lim inf — log uj 

£->■ 0 + ra->-oo tn 


X € X 


asix)ds G Js,e 


JO 


= lim liminf — loga;(Byi(s,5)) 

£ —^ 0 '(' ^ ^00 tfi 

> -Hs), 


where s G K is any point. A well-known (and simple) argument implies the required lower bound ( 68 ) for 
any interval C M. Thus, we need to establish (69). To simplify the notation, we shall consider only the 
case when s > s+ (assuming that s+ < oo). 

Step 2: Shifted measures. Let us fix s > s+ and denote et(a) = et{—a) and e{a) = e{—a). Since is a 
monotone increasing function mapping the interval —Jt„ = (—1 — onto (— cxd, oo) (see (46)), for 

any n > 1 there is a unique number G — such that ej (a„) = t„s. Following a well-known idea in 
the theory of large deviations, let us dehne a sequence of measures i/n on X by their densities 

= exp(-a„£^^^|^ - et„(a„)). 


Suppose we have proved that 

liminf e)) > 0. 

n—>-oo ' ' 

In this case, assuming that Q!„ > 0, we can write 


oj{Bn{s,e)) = / exp{aniujt ju + et„(an))di^n 

J Bn{s,£) 

> exp(f„a„(-.s - e) -f et„(an))r'„(B„(s,e)), 


(71) 


whence it follows that 


liminf — loga;(S„(s, e)) 

n—>-oo tfi 


> liminf 

n—foo 



e) + 


I 

tn 


(72) 


If we know that 


lim an = 

n—foo 


liminf — et„(a„) > e+, 

n—>-oo tn 


(73) 


then Q!„ > 0 for n large enough and inequality (72) and relation (70) immediately imply the required 
result (69). Thus, we need to prove (71) and (73). 


^In the formulation of Theorem 4.65 in [.lOPP], it is required that the limit of f „ ^ et^ (a) as n —> oo should exist for any a in the 
closure of J. However, the same proof works also in the case when the limits exist only for a a J. 
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Step 3: Proof of (73). Since S and 5t^ —^ 5, the first relation in (73) will be established if we 
show that 

liminfa„ = (5. (74) 

n—foo 

Suppose this is not the case. Then there is e > 0 and a sequence +oo such that — 1 < < 5 — e, 

where the first inequality follows from the fact that (a„) > 0 and e't^(-l) < 0. To simplify notation, 
we assume that the entire sequence {a„} satisfies this inequality. It follows that 

s'*" < s = —(«„) < —((5 — e) for any n > 1. (75) 

Since :^et^ (a) are convex functions converging to the smooth function e{a) for a G —J,hy Theorem 25.7 
in [Rock], we have 

lim — e[ {a) = e'{a) foranyae—J, 

n—>00 tn "■ 

and the limit is uniform on any compact subset of — J. Comparing this with (75), we see that s+ < e'{5—e). 
It follows that e' is constant on the interval — e, j] and, hence, by analyticity and the first relation in (15), 
the function e(a) vanishes. This contradicts the assumption that uj+{a) > 0 and proves (74). 

We now establish the second relation in (73). For any 7 S (0, S), we have 

et„ (a„) = et„ ( 7 ) + / (a) da > et„ ( 7 ) + (a„ - 7 ) 6 *^ ( 0 ), 

where we used the facts that e' is nondecreasing and that an > J for sufficiently large n > 1 , in view of 
the first relation in (73). It follows that 

lim inf 7 -et„(a„) > 6 ( 7 ) + ((5 - 7 )e'( 0 ). 

n—>00 tn 

Passing to the limit as 7 —5, we obtain the required inequality. 

Step 4: Proof of (71). Let us introduce trace class operators 

Qn=D^I^Tt^D^I^, Mn=tf\l-anQn)-^Qn, U > 1. 

Since a„ € —Jt„^ the operator I — anQn is strictly positive and, hence, invertible, so that Mn is well 
defined. Suppose we have shown that 

l'n{f{X„))= p,{f{Yn)), Xn = ¥„ = , MnX), U > I, (76) 

where / : K —K is an arbitrary bounded measurable function and p. is the centered Gaussian measure 
on X with the covariance operator /. In this case, taking / to be the indicator function of the interval J7s,e, 
we can write 


Vn{Bn[s,£)) = p{{x G X| \Yn{x) - s| < e}) =: Pn(£) 

Thus, the required assertion will be established if we prove that 

inf p„(e) > 0 for any e > 0. 

n>l 


for any n > 1. 


(77) 


To this end, let us assume that we have proved that 

Wl := sup ||M„||i < oo, tr(M„) = 2s. (78) 

n>l 


We now use the following lemma, whose proof is given in the end of this subsection (cf. Lemma 2 in [BD, 
Section 3].) 
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Lemma 4.2 Let /i be the centered Gaussian measure on X with the covariance operator I. Then for any 
positive numbers k and e there is p(k, e) > 0 such that 

p.{{x S X I |(a;, Mx) — tr(M)| < e}) > p{k, e) (79) 

for any selfadjoint operator M G f' satisfying the inequality |i-M|| i < n. 

In view of (78), we have 

Yn{x) - S = [x, \MnX) - 

Applying Lemma 4.2 with k = 2 we see that (77) holds. Thus, to complete the proof of the theorem, it 
remains to establish (76) and (78). 

Step 5: Proof of the auxiliary assertions. Simple approximation and analyticity arguments show that, to 
prove (76), is suffices to consider the case in which f{x) = where 7 € K is sufficiently small. Thus, 
we need to check that 

=p.{e^^^). (80) 

Recalling the construction of a„ and using the relation et (a) = — ^ log det(/ — aQt) (see (46)), we write 

= / exp(-(7f;;;^+a„)C^^|^ - e(„(a„))a;(dx) 

Jx 

= exp(et„( 7 f“^ + a„) - et„(Q;„)) =det(/- 7 M„) 

This expression coincides with the right-hand side of (80). 

Finally, to prove (78), we first note that the equality follows immediately from the choice of and the 
relation e^{a) = — aQt)~^Qt) ■ To establish the inequality, we start by using (39) and (63) to get 

the bound 

i-tn n/f2 

llQnlll < / < -fniklll. (81) 

Jo m 

Writing the spectral decomposition of the compact self-adjoint operator M„, we easily show that 


Y CtnQn ) Qn ! 

where A+ and A~ stand the positive and negative parts of a selfadjoint operator A, and we used that fact 
that > 0 for sufficiently large n (see (74)). Combining this relation with (81), we derive 

tr(M-) = tfhidl + anQf)~^Qf) < —Iklli- 
Recalling the second relation in (78), we conclude that 

f M'^ \ 

||M„||i = tr(|M„|) = tr(M„ -f 2M„ ) < 2 ( s H-||<^||i I . 

\ m J 

The proof of Theorem 2.5 is complete. □ 

Proof of Lemma 4.2. We set Y (x) = (x, Mx) and note that p,(Y) = tr(M). Let us denote by {P/, / C K} 
the family of spectral projections for M and, given a number 6 > 0, write M — M-^ + where 

= MP[_e, 9 ]. Accordingly, we represent Y in the form 

Y{x)=Y^^{x) + Y>\x), Y^\x) = {x,M^^x) 

Now note that the random variables Y-^ and are independent under the law p. It follows that the 
probability P{M, e) given by the left-hand side of (79) satisfies the inequality 

P(M,e) > p{{\Y>'^\ < e/2, < e/2}) = p{{\Y>^\ < e/2})p{{\Y^^\ < e/2}). (82) 
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We claim that both factors on the right-hand side of this inequality are separated from zero. Indeed, to 
estimate the first factor, we note that 

K > ||M||i > 6lrank(M>®), (83) 

where rank(M>®) =: Ng stands for the rank of Denoting by Xj the eigenvalues of M indexed in 

the non-increasing order of their absolute values, we see that 


|y>^(x))| = 


Ng 


i=i 


Ng 

< 

i=i 


- 1 |, 


where {xj} are the coordinates of x in the orthonormal basis formed of the eigenvectors of M. Combining 
this with (83), we derive 


^ Ne Ne 

^l{\Y>\x))\ < e/ 2 } > -^\ < y\ - < i^KNe)-^e} > 

^^ ' i=i 

where S = e 6 */( 2 K^), and p{S) > 0 is the probability of the event \x‘^ — 1 | < (5 under the one-dimensional 
standard normal law. To estimate the second factor in (82), we use the Chebyshev inequality; 

^l{\Y^\x)\ < e/ 2 } = 1 - p{Y^<>{x) > e/ 2 } - p{-Y^<^{x) > e/ 2 } 

> 1 - Ai(exp( 7 y^® - 7 e/ 2 )) + ^(exp(- 7 r^® - je/2 )), (84) 

where 7 > 0 is sufficiently small and will be chosen later. We have 

p{exp{'yY^^)) =exp{- 7 tr(M^^®) - i logdet(j - } 

= exp{-itr( 27 M^® -f log(/- 27 M^®))}. (85) 

Now note that if 4 | 7|0 < 1, then 

+ log(/ - 2jM^^) = ^ ^ ■ 

n—2 

Recalling that ||M-®|| < 9 and ||M-®||i < k and using the inequality |tr(ylB)| < ||A||i||i3|j, it follows 
that 

00 

|tr( 27 M^® -f log(/ - 2-fM^^)) \ < ^ | 276 »r-i 2 | 7 |/c < 

n—2 

Substituting this into (85), we see that, if I 7 I < (40)“^, then /i(exp( 7 y-®)) < exp( 4 K 7 ^ 0 ). A similar 
estimate holds for ^(exp(— 7 y-®)). Combining these inequalities with (84) and choosing 7 = we 
derive 

p{\Y^^{x)\ < e/2} > l-2exp{4K'y^9-'ye/2) = 1 - 2exp(-g|^). 

The right-hand side of this inequality can be made greater than zero by choosing a sufficiently small 0 > 0 
which will depend only on k and e. □ 


4.7 Proof of Theorem 2.6 


The proof of this result is verty similar to that of Theorems 2.4 and 2.5, and we shall only outline the proof. 
Part (1) Follows from Holder’s inequality as in the proof of Proposition 2.3 (2). 

Part (2) Since 0 S J/”, the fact that j/" is an interval follows immediately from the following property: if 
a G j/~, then 9a G J/" for 9 G (0,1). To prove the analyticity, note that, by Eq. (45), one has 


— Oit. 




(det(/ + f2rt)) ^ 

det(J — aD+Tt) 


• dw 
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This relation implies that the function 

et+(a) = logdet(/ + D+Tt) - ^ logdet(/ - aD+Tt) (86) 

= logdet(/ + i logdet(/ - 

1/2 1/2 

is real analytic in a on the open interval defined by the condition I — aD_^ > 0 and takes the 

value +00 on its complement (where the intersection of the spectrum of I—aD]^'^TtD]^^ with the negative 
half-line is nonempty). The above inequality coincides with the one defining . 

Part (3) The fact that is an interval follows immediately from its definition. To prove that j/" 
note that, in view of Hypothesis (G2), for any f, a S K we have 


I- 


aD]/^TtD 


1/2 

+ 






|q| 

6 + 1 ■ 


This expression is positive for | q ;| <6. 

To prove the existence of limit (27) and its analyticity on J^, we repeat the argument used in the proof of 
Theorem 2.5 (2). Namely, let us introduce the family of operators (a) = — aTt)~^, which are 

well defined for a € (—6, 6). Then the following analogue of relation (65) is valid; 


^D+(a)if) 




for / e L\X, da;^+(„)). 


The argument used in the derivation of (66) gives that 


^et-r(a) 



J 0 


while Hypothesis (G2) and the relation = £>+ valid for r G K imply that 


^—tsC 


Dt 


(7)e 


— tsC* 


= {D-^^ - - DZl)y" <m(i- 


Following again the argument in the proof of Theorem 2.5 (2), for a G (—6,6) we derive 

1 /■“ 

e+(a) = lim -et+(a) = - Jcr) d 7 . (87) 

t->-oo t Jo 

Now note that = DD^-d, whence it follows (cr) = ojjj (cr o d) = —ujjj (cr). Substituting this 

into (87) and recalling (67), we see that 

PCX 

e+(a) = / uijy {a)d'y = e{a) for a G (—(5, (5). (88) 

Jo 

We have thus established the existence of limit (27) on the interval {—5,5) C The fact that it exists for 
any a G J~^ and defines a real-analytic function can be proved with the help of Vitali’s theorem (cf. proof 
of Part (2) of Theorem 2.5). Finally, relation (28) is established by the same argument as (23). 

Parts (4-6) The proofs of the large deviation principle, central limit theorem, and strong law of large 
numbers for the time average of the entropy production functional under the limiting law a;+ are exactly 
the same as for lu (see Parts (3-5) of Theorem 2.4), and therefore we will omit them. 

Parts (7) The fact that the functions e+(a) and e(a) coincides on the intersection J’^’n J follows from (88) 
and their analyticity. The equality of the corresponding rate functions on a small interval around u;+(cr) is 
a straightforward consequence of (88) and the definition of the Legendre transform. □ 
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